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Motivation  for  this  Work 


We  are  interested  in  constructing  numerical  methods  for  constructing  non-hydrostatic 
mesoscale  and  global  atmospheric  models  (for  NWP  applications);  this  is  a  unified 
model.  The  reason  for  this  is  economics  -  one  (production)  model  is  cheaper  to  support. 

Currently,  in  the  U.S.  there  is  a  movement  to  construct  one  NWP  model  (NWS,  Navy, 
and  Air  Force).  This  National  Board  (NUOPC=National  Unified  Operational 
Prediction  Capability)  aims  to  develop  a  new  model  that  is: 

1 .  Highly  scalable  on  current  and  future  computer  architectures 

2.  Global  model  that  is  valid  at  the  meso-scale  (i.e.,  non-hydrostatic) 

3.  Applicable  to  medium-range  NWP 

4.  Applicable  to  decadal  time-scales 

The  following  talk  outlines  a  model  development  effort  to  meet  these  needs. . . 


Talk  Summary 


•  Governing  Equations 

•  Spatial  Discretization 

•  Preliminary  (Validation)  Results 

•  Parallel  Implementation 

•  Closing  Remarks 


Governing  Equations 

(compressible  Euler  equations) 
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Spatial  Discretization 


Primitive  Equations: 


dq 

dt 


+  V-F  =  S(q) 


mn 

Approximate  the  solution  as:  qN=^'l',<h  Fv  =  F(V/v)  SN=S(qN) 

i= 1 


Interpolation  O(N) 


Write  Primitive  Equations  as: 


R(qN)^  +  V-¥N-SN=£ 


Weak  Problem  Statement:  Find 
qN  e  X(Q)  V  y/  e  S 


2  =  {vreF/1(D):^ePA,(De)VQe}  (CG) 
Z  =  {y/eL2(Q.):y/ePN(ne)VQe}  (DG) 


such  that 

•  Integration  0(2N) 
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Spatial  Discretization 
(Comparison  of  CG/DG  Methods) 


Continuous  Galerkin  Methods 


Discontinuous  Galerkin  Methods 


High  order  accurate  yet  local  construction  (via 
DSS) 

Simple  to  construct  efficient  semi-implicit 
time-integrators 

In  high-order  mode,  primarily  used  with  quads 
and  inexact  integration  (e.g.,  using  Lobatto 
points  avoids  non-diagonal  mass  matrix  with 
slight  error  since  integration  is  0(2N-1)) 

No  analog  of  Lobatto  points  exist  on  the 
triangle  so  costly  to  use 

Excellent  scalability  on  MPP 


High  order  accurate  and  completely  local  in  nature 
(no  DSS  required  as  in  CG) 

High  order  generalization  of  the  FV  (but  with 
compact  support) 

Upwinding  and  BCs  implemented  naturally  (via 
Riemann  solvers) 

Not  so  easy  to  construct  efficient  semi-implicit 
time-integrators,  due  to  the  difficulty  in  extracting 
the  Schur  complement 

Since  matrices  are  all  local,  using  quads  or 
triangles  is  straightforward  and  one  need  not  worry 
as  much  about  exact  vs.  inexact  integration 

Excellent  scalability  on  MPP 


Preliminary  Results 
(Model  Description) 

•  Basis  functions:  3D  tensor  products  of  Lobatto-Gauss-Legendre  (LGL)  points. 
Elements  are  hexahedra  (Triangular  prisms  coming  soon). 

•  Time-Integrators  are:  explicit  SSP-RK,  IMEX-BDF2  (Schur  and  No  Schur),  Fully- 
Implicit  BDF2  (JFNK),  IMEX-RK  (currently,  No  Schur  only) 

•  Mesoscale  (limited  area)  and  Global  (spherical  domain)  options 


Mesoscale 


Global 


Preliminary  Results 

(Linear  Hydrostatic  Ridge  and  Mountain) 


•  Flow  ofU=20  m/s  in  an  isothermal  atmosphere. 

•  LH  Ridge:  Witch  of  Agnesi  ridge:  Mountain  height  =  1  m  with  radius  10  km. 

•  LH  Mountain:  Solid  of  revolution  of  Witch  of  Agnesi:  Mountain  height  =  1  m  with 
radius  10  km. 

•  Absorbing  (sponge)  boundary  condition  implemented  on  lateral  and  top  boundaries. 


LH  Ridge 


LH  Isolated  Mountain 


Linear  Hydrostatic  Ridge 
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Preliminary  Scaling  Experiments 
(Performed  on  Ranger  TACC) 


32x32x32  elements  with  4th  Order  48x48x48  elements  with  4th  Order 

Polynomials  Polynomials 

(2  Million  Grid  Points)  (7  Million  Grid  Points) 


A  Multitude  of  Challenges  Remain 


•  Further  dry  physics  validation  is  necessary  (e.g.,  Baroclinic  Instability  problems). 

•  Simple  moisture  has  been  tested  in  2D  (manuscript  almost  finished)  and  now 
implementing  it  in  3D. 

•  Full  sub-grid  scale  parameterization  needs  to  be  included  (can  compare  against 
older  hydrostatic  version  called  NSEAM). 

•  Interesting  question  is:  how  will  the  NH  and  H  models  compare  in  terms  of  both 
solution  quality  and  cost? 

•  Adaptivity  will,  eventually,  be  included  (as  in  A.  Muller)  but  I  envision  only 
using  triangular  prisms. 

•  Explicit  scalability  is  great  but  must  improve  on  Semi-Implicit  performance 
(different  time-integrators  and  new  approaches  for  DG  such  as  in  M.  Restellf  s 
talk  on  hybridized  DG). 


